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We consider tlie electronic structure of, and magnetic exchange (spin) interactions between, nom- 
inally nonmagnetic Eu^"*" ions (4/^^, S=3, L=3, J=0) within the context of the rocksalt structure 
compounds EuN and EuP. Both compounds are ionic [Eu''"'"; A''''" and P^~\ semimetals similar to 
isovalent GdN. Treating the spin polarization within the 4/ shell, and then averaging consistent 
with the J=0 configuration, we estimate semimetallic band overlaps (Eu 5d with pnictide 2p or 3p) 
of ~ 0.1 eV (EuN) and ~ 1.0 eV (EuP) that increase (become more metallic) with pressure. The 
calculated bulk modulus is 130 (86) GPa for EuN (EuP). Exchange (spin-spin) coupling calculated 
from correlated band theory is small and ferromagnetic in sign for EuN, increasing in magnitude 
with pressure. Conversely, the exchange coupling is antiferromagnetic in sign for EuP and is larger 
in magnitude, but decreases with compression. Study of a two-site model with Si -82 coupling within 
the J = 0, 1 spaces of each ion illustrates the dependence of the magnetic correlation functions on 
the model parameters, and indicates that the spin coupling is sufficient to alter the Van Vleck sus- 
ceptibility. We outline a scenario of a spin-correlation transition in a lattice of5' = 3,L = 3,J = 
nonmagnetic ions. 
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I. INTRODUCTION 

The readily observable behavior of the angular momen- 
tum and associated magnetic moment of rare earth ions 
is one of the more obvious successes of the quantum me- 
chanical description of atomic magnetism. Hund's rules 
for the total angular momentum J = L + S in terms of 
orbital filling (giving orbital L and spin S angular mo- 
menta) are simple and successful. The most striking pre- 
diction is when the angular momenta in a spin polarized 
S and orbitally polarized L ^ ion conspire to give 
a nonmagnetic ground state J = 0. The canonical exam- 
ple is Eu^+f: 5* = 3,L = 3, J = 0. Whereas Eu^+ is 
nonmagnetic experimentally (unresponsive to magnetic 
fields, treatment of the high spin polarization S*^ — 12 
(we use ?i=l) is essential for obtaining any reasonable 
model of the electronic structure. 

The orbital polarization and associated moment rfiL — 
HbL is a secondary effect as far as the determination of 
the electronic structure is concerned, because the mag- 
netic coupling proceeds via the effect of the spin moment 
rhs = 'ifj.BS on the electronic structure and on neigh- 
boring spins (via exchange coupling). L is however cen- 
tral in determining J and the total moment, which for 
Eu3+ is large: \M\ = ^ib\{L + 25^)1 = mb|(J -f ^)| = 
^ib\S\ ~ 3.5fiB i{M^)j=o = 12/i|)- Of course, in the 
J=0 ground state the moment has zero projection along 
any axis including any applied field, hence the "nonmag- 
netic" character of the ion. These considerations do not 
change the fact that the ion has a large intrinsic spin 
polarization, and that there will be exchange (spin-spin) 
coupling between neighboring spins; {Si)=0 does not im- 
ply {Si ■ Sj)=0. Such couphng may be difficult to ob- 
serve because the most direct means of observing mag- 



netic coupling is via response to a magnetic field. 

In this paper we address some of the questions raised 
by the interaction between these strongly polarized yet 
nonmagnetic ions within the context of two of the sim- 
plest solid state realizations, the binary compounds EuN 
and EuP. EuN is a semimetal with a simple rocksalt 
structure. Although it has been known for fifty years, 
there is little characterization of this compound in the 
literature. Reported values of the lattice constant are 
consistent at 5.020±0.006 AiSi^. Only a very small 
non-stoichiometry range exists, although magnetic Eu^+ 
impurities can occur in sufficient concentration in some 
samples to mask the susceptibility of the Eu3-f ion in 
EuN.^ There is also sparse literature for EuP, a com- 
pound both isostructural and isovalent to EuN. The lat- 
tice constant of EuP is measured^ to be 5.756 A. The 
15% increase in lattice constant (50% larger volume) 
compared to EuN has a significant effect on the exchange 
interaction, probably by shifting the balance between 
competing mechanisms of exchange coupling. 

Although any projection (Si) vanishes, when two 
neighboring ions are coupled the spin-spin correlation 
{Si ■ S2) should be nonzero. Correspondingly there will 
be a term in the microscopic Hamiltonian K X]<y> ' 
for all interacting pairs. (We denote the exchange cou- 
pling by K and retain J for other uses.) This coupling 
in itself breaks the spherical symmetry of the ion, cou- 
pling the J = 0, 1, 2,.., L+S states and allowing mag- 
netic behavior to emerge. This effect is already seen 
in the magnetic susceptibility of the Eu'^^ ion, where 
the magnetic field breaks the spherical symmetry, mixes 
J—1 with J—0 character, and gives rise to the van Vleck 
susceptibilitji£i2i2ii£. One may consider then whether this 
exchange coupling results in a phase transition; that is, 
a coupling strength or temperature for which the spin 
correlation length diverges. If the J = state of each ion 
is predominant, there will be only a weak magnetic sig- 



nal and also minor entropy [A:Blog(2J+l)^0]; this is the 
observed nonmagnetic behavior of the Eu'^+ solid. Van- 
ishing entropy implies that ordering could not be seen 
thermodynamically; on the other hand, with free energy 
devoid of magnetic entropy ordering could occur at much 
higher temperature. Neutron scattering, which has sensi- 
tivity to the difference between spin and orbital moments, 
should be able to detect long range anti-aligned ordering 
of the spins (hence of the magnetic moments) because of 
the AFM Bragg peak. This still primarily J=0 lattice of 
ions may therefore show a new kind of phase transition 
in which the {Si ■ Sj) correlation length diverges but with 
small (perhaps negligible) decrease in entropy. As the 
J=l state gets mixed in, by exchange coupling and by 
crystalline anisotropy, magnetic behavior will assert itself 
and the transition would become easier to detect. This 
process can be controlled, and the tendencies enhanced, 
by pressure-induced volume reduction. 

In this paper we begin to pursue this line of reasoning 
quantitatively. We first apply density functional based 
correlated band theory (LDA-I-U) to study the electronic 
structures of EuN and EuP and their equations of state. 
In the process we investigate the extent to which corre- 
lated LDA-I-U results can reproduce Hund's rules (more 
specifically, the z-components of the angular momenta) . 
We then apply conventional procedures to express the 
coupling in terms of a Heisenberg exchange coupling 
constant K between nearest neighbors. Using this ex- 
change coupling, the known spin-orbit coupling (SOC) 
constant, and a single-ion anisotropy parameter, we look 
at a toy model of coupled ions, each of which has the 
S = 3,L = 3,J — ground state of Eu'^+. Specifically, 
we consider whether spin-spin correlations between these 
J=0 ions can be detected, and whether they may in fact 
lead to a phase transition (divergent correlation length). 



II. CALCULATIONAL METHODS 

For Eu as for most rare-earth ions it is necessary to 
use the LDA+U method (or another orbital-functional 
approach) to separate the occupied 4/ orbitals from the 
unoccupied ones, because LDA invariably results incor- 
rectly in narrow partially filled 4/ bands straddling the 
Fermi level. For our LDA-I-U calculations, we have cho- 
sen a value U= 8 eV for the effective Coulomb repulsion 
and a value of eV for the intra-atomic 4/ exchange 

interaction on the Eu ion. These values were chosen by 
comparison with other Eu compounds. ^^'^^ Since U is 
so large, the specific value of U should not substantially 
change the conclusions about magnetic coupling which 
are obtained. The value of J''^ is typical of what is used 
for rare earth ions, and the effects of changing it have 
been found to be small. The non-intuitive effect on the 
electronic structure of the "exchange" as it arises in 
the LDA+U method will be discussed in a Sec. Ill C. 

All band structure calculations were performed using 
the program^. Wien2k, which employs an APW-l-lo (aug- 
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FIG. 1: The AFMII spin ordering of a rocksalt structure com- 
pound, such as is common in the transition metal monoxides. 
All second nearest neighbor spins, which are coupled by the 
180° superexchange coupling, are anti-aligned in this ordered 
phase. 



mented plane wave plus local orbitals) basis set. This 
allows for a lower plane wave cutoff than is necessary 
for LAPW calculations, and consequently the KK^ax 
was set to 7.00. We used the LDA exchange-correlation 
parameterization of Perdew and Wang^- and added a 
Hubbard U according to the prescription of Anisimov 
et aiiS*i&. In this scheme, the double counting terms 
are subtracted assuming a fully-localized (i.e. atomic) 
limitpii The APW sphere values were set to 2.1 a.u. for 
Eu and 1.8 a.u. for N, in EuN and to 2.4 a.u. for Eu and 
2.2 a.u. for P in EuP. These choices resulted in negligible 
core charge leakage in both cases. 

Band structure methods are seriously limited in the 
treatment of rare earth ions due to their classical treat- 
ment of the spin and orbital angular momentum (and for 
other reasons). As usual, we assign a direction (local z- 
axis) to 5, L and deal explicitly with only the correspond- 
ing projections Sz,Lz] one might say that we go into the 
reference frame of the spin on one Eu ion. The closest ap- 
proximation to the J—Q state is then Sz = 3, Lz = —3, 
which uniquely specifies the occupation of one-electron 
4/ orbitals. One challenge the LDA-I-U method faces is 
this: Hund's rules specify (1) coupling all spins to maxi- 
mum S, (2) coupling all orbital moments to maximum L, 
and only then (3) considering the effects of spin-orbit to 
obtain J. One electron approaches, on the other hand, 
typically incorporate spin-orbit coupling at the single 
electron level s + £ = j. 

Crystal field effects in both EuN and EuP are likely 
overestimated by the LDA, resulting in an orbital mo- 
ment that is "over-quenched" . Comparing energies of 
the Sz — S, Lz — —3 state and a mixed state with a di- 
minished orbital moment of L^ = -1.5, we find an energy 
difference of about 25 mRy for both EuN and for EuP. 
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The Jz=0 state is (meta)stable and at some volumes the 
system converges preferentiahy to this state, despite its 
somewhat higher total energy. Such energy differences in 
LDA+U sometimes seem to approximate excitations of 
the 4/ shell."- For these Eu pnictides and at the current 
level of understanding, these energy differences should 
probably not be taken too seriously; furthermore, the 
only noticeable changes in the band structure between 
the two orbital configurations occur in the splittings be- 
tween occupied 4/ states (see Section III C). We therefore 
use the lower energy state in our interpretations, keeping 
in mind that, in reality, the orbital moment is likely to 
be larger than our calculated moment. 

It is unclear from the outset whether the exchange 
coupling will favor alignment (ferromagnetism FM) or 
anti- alignment, AFMII, for example. (We use FM and 
AFM to refer to the relative spin alignment in our calcu- 
lations.) While the great majority of magnetic insulators 
are AFM, in particular the transition metal monoxides, 
the rocksalt structure EuO and Eu chalcogenides are FM 
or become FM under pressureiii In the calculations we 
present, we double the unit cell to investigate the ener- 
getics of AFMII ordering. To obtain an AFM solution, 
AFM symmetry of the wavefunctions and of the LDA-I-U 
occupation matrix has been imposed. Occasionally the 
AFM symmetry constraint was released to ensure that 
the solution remains stable without constraint. 



III. ELECTRONIC STRUCTURE 
A. EuN 

As mentioned above, LDA is qualitatively incorrect for 
the localized 4/ states of EuN, and LDA-I-U must be used 
to separate occupied from unoccupied 4/ states. In Fig. 
121 the LDA and LDA-I-U band structure results for the 
majority bands of FM EuN are displayed. The LDA re- 
sults are qualitatively incorrect and will not be discussed. 
The N 2p bands lie in the range -4 to eV. The lower con- 
duction bands are of Eu 5c? character, disrupted by the 
single empty majority spin / band, and overlap the N 2p 
bands at the X point of the zone to form a semimetallic 
band structure. The S=3 configuration is represented 
by six fiat 4/-bands in the -4.5 eV to -7 eV range sepa- 
rated roughly by U from the single unoccupied majority 
4/ state lying 1 eV above Ej?. The 4/ bands are so lo- 
calized that they do not overlap even with neighboring 
4/ orbitals, making direct exchange negligible. As we 
will show, the Eu 5(i states are considerably extended 
and provide one mechanism for coupling of the localized 
spins (and hence the moments) on the Eu ions. 

The effects of the spin polarization of the localized 4/^ 
shell on the electronic structure are considerable and non- 
intuitive. Via the on-site Hund's coupling (contained in 
the LDA exchange potential) the Eu hd bands are polar- 
ized with the same orientation, with an exchange split- 
ting of 0.6 eV. The N 2p bands are also polarized, but in 




FIG. 2: (color online) Band structures of EuN in the simple 
unit cell with the 4/ character highlighted by broadened sym- 
bols. Left panel: the incorrect metallic solution given by LDA. 
Right panel: The LDA-I-U solution has six filled majority 4/ 
bands and a single empty 4/ majority spin band, separated 
roughly by the value of U. The occupied 4/ band splittings 
(heavy lines) are discussed in the text. The Eu 5d conduction 
band dips below the Fermi energy to cross the valence N p 
bands at the X point; Eu 5d character is denoted by the (red) 
circles. 



the opposite direction. As a result, the majority 2p bands 
lie above the minority bands, and the induced magneti- 
zation on N is negative. While the majority Eu bd band 
overlaps the majority 2p band at Ei?, there is a gap be- 
tween the respective minority bands. Hence FM EuN is 
half metallic within conventional LDA-I-U band theory. 

A comparison of the FM band structure of EuN in a 
compressed face-centered cubic (FCC) unit cell and the 
AFM band structure in a similarly compressed, doubled, 
rhombohedral unit cell is shown in Fig. |3| Compression 
moves the Eu 5d bands down relative to the N 2p bands, 
increasing the filling of bd states at the X point {L in 
the rhombohedral cell) and beginning to fill a minority 
Eu bd pocket. Thus half metallicity is lost under moder- 
ate compression (pressure). The compensating holes go 
into the 2p band, at F and at the X point. The remain- 
ing complete spin polarization of the N 2p holes may be 
important for the 4/ — 4/ exchange processes (Sec. IV). 

The only previous calculations on EuN were by 
Horne et ali^ using the self-interaction-corrected LDA 
(SIC-LDA) approach implemented within the linearized 
muffin-tin orbital band structure code. Horne et al. re- 
ported only the DOS for FM spin alignment, with the 
results being significantly different from what we have 
obtained, as the considerable differences in the methods 
might suggest. Regarding the Eu bd and N 2p bands, 
they obtained a semiconducting result with a separation 
of 1.4 eV, whereas LDA-I-U gives an overlap of several 
tenths of an eV for FM alignment. The N 2p bandwidth is 
15% narrower than the LDA-fU value of 4.5 eV. IN SIC- 
LDA the unoccupied majority 4/ orbital lies at the very 
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FIG. 4: (color online) ie/t panel: bands of FM EuP in the 
FCC cell with majority spin character highlighted. Right 
panel: AFM EuP in the rhombohedral cell at equilibrium 
volume. Note that the 4/-3p mixing near Ef is so strong in 
the AFM case that the unoccupied 4/ band above Ef is not 
readily obvious. 



FIG. 3: A direct comparison of EuN FM and AFM band 
structures at a 4% compressed lattice constant (15 GPa pres- 
sure). The FM configuration (left) is shown in the FCC unit 
cell with all majority spin states highlighted, while the AFM 
configuration (right) is shown in doubled rhombohedral cell. 
The rhombohedral (shaded) Brillouin zone is shown inside the 
larger BZ of the FCC unit cell below. The important Eu d 
overlap with the N p bands occurs at the X point in the large 
BZ, corresponding to the L point of the small BZ. 



bottom this gap, mixes strongly with the N 2p states, ac- 
quires a 1 eV width and becomes slightly occupied. The 
SIC-LDA result thus is a half metal (for FM alignment). 
The minority 4/ bands (unoccupied) lie about 4 eV above 
Ei^, quite similar to what we find from LDA+U. The ma- 
jority 4/ states in SIC-LDA lie about 12 eV below Ep 
with 0.5 ev width, whereas those obtained from LDA-f U 
are centered only 5.5 eV below Ep with a 2 eV width. As 
we discuss in Sec. Ill C, the 4/ state splittings in LDA-f U 
are a result of the anisotropy of the exchange interaction. 
These various differences generally reflect the main dif- 
ference in the methods: SIC-LDA lowers (very strongly) 
the occupied 4/ states, whereas LDA-I-U lowers occupied 
and raises unoccupied 4/ states each by some fraction of 
U. 

Other rare earth nitrides have been reported in experi- 
mental studies to be semimetal o^°i^^i^^ as our results give 
for EuN, while the Group IIIA mononitrides (BN, AIN, 
GaN, InN) are known to be wide gap semiconductors. 
Calculations of the electronic structure of the closely re- 
lated compound GdN have been given by Petukhov et 
al.^ who used an atomic sphere approximation for the 
potential and treated the 4/ states corelike (thus removed 
from the band structure problem). These calculations 



gave a small (0.1 eV) overlap between the valence band 
maximum at T involving N 2p states, and a conduction 
band minimum at X comprised of strongly Eu 5c? charac- 
ter. This band structure is very much like what we obtain 
for EuN (as well as for antiferromagnetic EuP below) . By 
averaging over the spin splitting (and also surveying the 
AFM bands) we deduce a semimetallic band overlap of 
~ 0.1 eV for EuN, and - 1.0 eV for EuP, both of which 
increase under pressure. 



B. EuP 

The band structures of FM- and AFMII-ordered EuP 
are shown in Fig. ^ The main features of EuN remain 
in EuP, specifically the character of the valence and con- 
duction bands and the semimetallic band crossing at the 
X (L) point. In EuP, the valence P 3p bands that deter- 
mine Ep are higher in energy (relative to the 4/ bands) 
than the analogous 2p bands in EuN. Conversely, one 
may adopt the viewpoint that the Eu 5d and 4/ bands 
are shifted downward relative to the P 3p bands (see Fig. 
|3J| . One feature clearly relevant to magnetic interactions 
is 4/ — 3p hybridization. The unoccupied majority 4/ 
band sits just above Ep and mixes strongly with the P 
3p band near X along the T — X line. Holes in the 3p 
band are consequently near the F point only and not from 
the majority spin band near X which is pushed down by 
hybridization with the / band. In contrast to EuN, EuP 
with FM alignment is not a half-metal and significantly, 
the 3p holes are only weakly spin-polarized. 
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FIG. 5: (color online) A partial DOS plot for FM EuN and 
AFM EuP. Top panel: The orbitally resolved 5d, 4/ states of 
Eu and 2p states of N at a compression of 12% by volume 
(15 GPa). The unoccupied 4/ band is more than 1 eV above 
the filled N 2p bands. The bottom of the valence 5d band 
is filled, providing the mechanism for superexchange. Bottom 
panel: At equilibrium volume, the unoccupied / band sits im- 
mediately above Ef, and the occupied bands are completely 
separated from the filled P 3p bands. 



configuration L, S, J. Solovyev et al. have addressed this 
question;^ and Gotsis and Mazin have obtained encour- 
aging results^! for some Sm intermetallic compounds. 

The separations of the filled fiat 4/ bands in Figs. |21 
13 El are due to the strong orbital-dependence of the ex- 
change interaction, whose full anisotropy is incorporated 
into the LDA-I-U code. These differing energy shifts are 
determined by the Hund's exchange constant 3^. To 
clarify the origin of the splitting between the f-bands, we 
present an explicit example of the potential matrix ele- 
ments used in the LDA-I-U procedure for the case where 
the density matrix is diagonal, n'^^, = n'^Smm' (which 
avoids uninstructive complications that arise in the gen- 
eral case). Then the only matrix elements which enter 
the LDA-I-U energy are the direct and exchange integrals 
of the form 

Kmm' = {mm'\Vee\mm') (1) 
Jmm' = {mm' \Vee\m' m) 

In the atomic limit these are given simply in terms of the 
i = 3 Slater integrals Fq, F2, F4, Fg and angular factors. 
In the LDA+U method these are represented in terms of 
two constants: "direct" U and "exchange" 

U = F''- J" = (F2+^4)/14 (2) 
and the ratio ^4/^2 - 0.625 ^5/^2 - 0.494 



C. 4/ Occupation and Hund's 2nd rule 

It should be understood at the outset that any band 
picture of the 4/ states has serious intrinsic limitations 
that have been recognized for some time and occasionally 
addressed^ A band structure is hoped to represent, in 
a mean field manner, the electron removal spectrum for 
occupied states and the electron addition spectrum for 
unoccupied states. The open 4/ shell is a highly corre- 
lated local system: individual spins are strongly coupled 
and aligned to form the total spin 5, orbital occupations 
are correlated to give L, and Hund's third rule (SOC) 
gives the ground state J. In such ions there is a manifold 
of low energy excitations in which the particle number 
remains unchanged (J 7^ multiplets at higher energy). 
These multiplet excitations, often affected strongly by 
crystal fields, lie outside the realm of band structures, 
even in a 'zeroth order' description. 

Density functional theory (DFT), which gives rise in 
the process of energy minimization to the Kohn-Sham 
band structure, is itself a many body theory for the sys- 
tem energy. There have been a few indications that 
the LDA-I-U approach is useful in modeling (perhaps 
predicting) some of these particle-number conserving 
excitations»i2^ More specifically, one can ask whether 
the LDA-I-U approach is capable of reproducing Hund's 
rules, i.e. obtaining the lowest energy for the correct 



With U = 8 eV and J = 1 eV as we have used, the 
corresponding Slater integrals are: Fq = 8.0 eV, F2 = 
11.92 eV, F4 = 7.96 eV, and Fg = 5.89 eV. (Fq is strongly 
screened in the solid so atomic values cannot be used 
directly.) 

If is set to zero, these simplify to an isotropic 
repulsion lAmm' = Fq, and to a diagonal and orbital- 
independent J^mm' = Fq 5mm', the latter simply re- 
moving the self-interaction. Thus for J^=0 no distinc- 
tion between different orbitals is made, and filled orbitals 
will all be shifted downward equally, while unfilled or- 
bitals will be shifted upward equally. An important and 
unsatisfactory result is that the energy differences based 
on orbital quantum number m (which give rise to Hund's 
2nd rule splittings in the atomic limit) are absent when 
J»=0. 

With J-^ ^ the diagonal elements become or- 
bital dependent, and each occupied orbital obtains (self- 
consistently) its own downward shift in energy. To il- 
lustrate the strong orbital dependence (anisotropy) we 
present explicitly the Umm' and J^mm' matrices for Eu 
in Appendix I. The resulting differences in eigenenergies 
can be pronounced, as in Fig. ^Iwhere different occupied 
4/ bands are separated by as much as 2 eV. therefore 
has two roles: (1) it provides the anisotropy of the direct 
Coulomb repulsion Umm' (which is independent of m, m' 
if J^=0), and (2) it has the more commonly understood 
aspect of providing the intra-atomic exchange coupling 

\J m,m' 5 m, ^ ?7i . 
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D. Equation of State 

We fit the energy vs. volume curves for the FM/AFM 
spin configuration of EuN/EuP to an equation of state2&: 

E = a + hV-^'^ ^cV-^'"^ + cV-^ (3) 

and extracted the LDA equilibrium constants and bulk 
moduli. EuN has a calculated equilibrium volume of 
406.9 a.u.^ , 5% smaller than experiment (consistent with 
typical LDA error) and a bulk modulus, B=130 GPa. 
The volume of EuP is 612.6 a.u.^ , consistent with the 
larger P'^" ion, and again is 5% smaller than experiment. 
It has a bulk modulus of 86 GPa, making it much softer 
than EuN. The lattice constants of EuN and EuP have 
been calculated previously by Horne et al.^ whose SIC 
results are very similar to ours. 



IV. EXCHANGE ENERGY AND 
INTERACTIONS 

In the rocksalt structure each Eu ion has potential 
exchange coupling to its twelve nearest Eu neighbors 
through an intermediate pnictide ion at a 90° angle, and 
to each of its six second neighbors at a 180° angle. The 
Goodenough-Kanamori superexchange rules lead one to 
anticipate a ferromagnetic nearest neighbor interaction, 
while the second nearest neighbor interaction will be an- 
tiferromagnetic. These guidelines are based on pertur- 
bation expansions appropriate for insulators, and may 
be misleading here where other mechanisms also arise 
(RKKY interactions through the electron and hole car- 
riers, specifically). We probe some energetics of the su- 
perexchange processes by calculating the energy differ- 
ences between ferromagnetic and AFMII spin alignments, 
using the same rhombohedral cell to improve numerical 
precision. For the experimental lattice constant of EuN, 
there is a small energy difference Ai? — F,fm - ^afm 
= -4.9 meV. AE increases in magnitude with compres- 
sion and reaches a value of -10.1 meV when the lattice 
parameter is reduced by « 4%. In EuP, the energy dif- 
ferences are somewhat larger and of the opposite sign, 
with AE= 12.8 meV at equilibrium volume decreasing 
to AE = 9.2 meV at a compression by 4% of the lattice 
constant. The calculated ground state of EuN is there- 
fore FM, while that of EuP is antiferromagnetic. Our 
single total energy difference is only sufficient to give a 
single exchange constant. As can be seen in the equation 
below, the effective exchange constant (which includes 
all exchange interactions between all neighbors) may be 
attributed solely to the first neighbors or solely to the sec- 
ond neighbors, resulting in an identical expression. Us- 
ing nearest neighbors only in a classical S=3 Heisenberg 
model gives for the two magnetic states: 

Eafm = 0, Efm=QKS'\ AE = 6KS^. (4) 



Using second nearest neighbors only gives: 

Eafm = -5KS^, Efm ^ iKS"^, AE = QKS'^. 

Extracting the exchange constant from these formu- 
las gives KeuN — -0.09 meV (IK) at zero pressure and 
K_EtiAr = -0.19 meV (2.2K) at 4% reduced lattice constant 
or 15 GPa . For EuP, the exchange interaction is some- 
what larger, V^euP = 0.24 meV (2.8 K) at equilibrium 
pressure, and opposite in sign. It also has the oppo- 
site pressure dependency which, due to the smaller bulk 
modulus, produces a decrease of nearly 25% in KeuP at 
a pressure of 7.5 GPa (see Fig. |HJ|. We observe that 
for the larger compound, EuP, the effective interaction 
between ions decreases with increasing pressure despite 
the decreased distance between neighbors, possibly due 
to decreasing competition from first neighbor exchange. 
EuN is firmly within the FM regime and pressure has the 
expected effect of increasing the coupling. Recent pres- 
sure studies show that EuN retains the rocksalt structure 
up to 70 GPaiS Assuming that the lattice constant and 
K maintain a linear relationship at high pressure (Fig. |5| 
indicates that in reality it is slightly sub- linear), we can 
extrapolate that K = 0.3 meV (3.5 K) at 70 GPa. As 
we will show in the next section, an exchange coupling of 
this magnitude can have observable effects on the suscep- 
tibility. Despite a difference in Eu valency, these results 
arc quite consistent with the trends found for the series 
of Eu^+X compoundsii, namely that compounds with a 
large lattice constant are AFM, while compounds under 
high pressure or with smaller lattice constants are FM. 
The exchange constants of Ref. [l^ were resolved into 
Ki and K2 separately, both of which are comparable in 
magnitude to our single effective K. 

V. A TOY MODEL STUDY 

Since the states of ions are most readily described 
in terms of their total angular momentum Jj it is nat- 
ural and customary to use the effective Hamiltonian 
KeffJi ■ J2 for the coupling. However, the actual ex- 
change interaction proceeds via the spins Si and 6*2, so 
we begin this investigation by considering the more mi- 
croscopic model Hamiltonian expressed in terms of Si 
and Li. Our two identical site Hamiltonian (i=l, 2) is 

i/ = A-Si-S2+A^L,-S,-^M,-B+^(D-L,)^ (5) 

i i i 

This Hamiltonian includes the interactions we expect 
to be important in the investigation of magnetic behav- 
ior in EuN or EuP (although interactions may extend 
to more neighbors). We choose the spin-orbit constant 
to be A=323 cm^^ — 40 meV, a typical value for the 
''Fq -^'^ Fi transition observed by spectroscopioSiSSiSiiS 
measurements on several Eu^^ compounds. The ex- 
change parameter K will be varied to examine the macro- 
scopic effects of the exchange coupling, with the physi- 
cally relevant values taken from our LDA+U calculatios. 
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FIG. 6: The inferred exchange constant, K, for EuN and EuP. 
top panel: K for EuP is positive (antiferromagnetic order) 
and decreases with pressure. In analogy with the 3d transi- 
tion metal monoxides, AFMII order may more likely be due 
to next-neighbor AFM coupling than to near-neighbor AFM 
coupling, lower panel: K for EuN is negative (ferromagnetic 
order) and increases in magnitude with pressure. 
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FIG. 7: Alignment of spins and magnetic moment as a func- 
tion of K/\. The horizontal line in each case marks the value 
of maximum possible alignment possible in the restricted 
Hilbert space of the model 



The most conservative estimate is obtained by directly 
importing the K derived by mapping energy differences 
onto the Heisenberg model. This gives an upper bound 
of K/\ ^ 0.01 for EuN (under high pressure) and for 
EuP at equihbrium volume. It may be more realistic to 
establish K by directly forcing the toy model to repro- 
duce the density functional AFM-FM energy differences 
(per Eu ion pair). This gives a much higher value ratio 
of K/\: 0.21 for EuP and 0.08 for EuN, both at their re- 
spective equilibrium volumes. The single-ion anisotropy 
D is taken to be of similar magnitude (2 meV), with the 
easy axis defining the z axis. 

Our basis includes J=0 and three J=l states on each 
of the two sites. The next J multiplct, ''F2, is 2A ~ 900A' 
higher in energy than ''i^i, so it will not contribute to the 
low energy regime we explore. This basis of four states 
per site leads to a 16 x 16 Hamiltonian matrix. Matrix 
elements are readily evaluated in the \L, L^', S, Sz) basis; 
the expansion of the states | J, Jz) states {J — Ji + J2) in 
the \L, Lz] S, Sz) states is given in Appendix II. 

For K=0, D=0, Hund's rules give the Ji—0 ground 
state on each site (and trivially Jtot—O)- As discussed 
in the Introduction, any projection of each magnetic 
moment Mi = fisiLi + 25*^), and therefore the total 
M = Ml + M2, vanishes for uncoupled spherical ions 
although its mean square value is large. Exchange cou- 
pling K introduces mixing between between the states of 
the two sites, but preserves the spin rotational symmetry 
of the Hamiltonian. Expectation values of both individ- 
ual and total magnetic moments (and individual angular 
momenta) therefore are zero, but correlation functions 



such as {Si ■ S2) and (A/i • M2) grow as the strength of 
the coupling grows. Single-site anisotropy D 7^ or field 
B destroys the rotational symmetry of H and introduces 
additional mixing of states. 

Fig. [7| shows the increase of the spin-spin and mag- 
netic moment correlation functions in the ground state 
as K increases. In the limit K/X —^ 00, the model reduces 
to a two site Heisenberg model, with a well known ana- 
lytic solution: Si ■ S2 = {S^ - Sf - S'|)/2 = -12, where 
5 = in the ground state and Si—Li=3. The restricted 
basis set we use does not span the entire Hilbert space 
of the Heisenberg Hamiltonian and the spin-spin correla- 
tion function therefore asymptotes to a smaller value of 
-7.18. The higher J multiplets necessary for (^i • S2) to 
attain the analytical ground state value are moved too 
high in energy by the scale of the spin-orbit parameter, 
therefore -7.18 should be considered the maximum anti- 
alignment of spins within this system. The asymptotic 
and analytically calculated expressions for (Mi • M2) are 
also different because of the restricted basis set. 

For small K the order parameter grows linearly with 
K with steep slope, and at K/\ — 0.175, {Si ■ S2) reaches 
half its maximum value, as shown in Fig. [7| A simi- 
lar conclusion holds for the moment-moment correlation 
function. Rotational symmetry of the system is bro- 
ken, allowing a non-zero moment and ready detection, 
by single-ion anisotropy D ^ or by application of a 
magnetic field B. 

A calculation of vs. B in thermal equilibrium, 

presented in Fig. |H1 shows two clear crossovers at fields 
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FIG. 8: Magnetic moment vs. applied magnetic field for a 
thermodynamic ensemble of two site "molecules" described by 
Eq. 1^ The curves are sensitive to the value of the exchange 
constant K, suggesting that this quantity can be probed with 
experimental measurements of M vs. B. 



of ^ 720 and 1450 Tesla. These transitions exist only for 
relatively high values of K and low temperatures and cor- 
respond to qualitative changes in eigenvectors as the ex- 
ternal B-field overcomes the exchange and spin-orbit en- 
ergies respectively. At higher temperatures or lower K/A, 
every state in the Hilbert space is present even in the 
ground state and the curve evolves smoothly. Because 
fields of this magnitude are experimentally unattainable 
and would be strong enough to mix in even higher J mul- 
tiplets (J=2,3...), only the low field part of the curve in 
Fig. IHlcan be taken seriously. In low fields, (M) is linear 
in B with a slope (susceptibility x) that is strongly de- 
pendent on the exchange coupling. This sensitivity sug- 
gests that the degree of spin-spin correlation in a system 
containing Eu'^"'" ions can be probed by examining the 
temperature dependence of the magnetic susceptibility. 

Fig. El shows the calculated x(T; K, A) curves for 
different values of K within our two site model. For 
small K, the curves show typical Van Vleck paramag- 
netic behavior—, characterized by an initial flat region 
which quickly decreases towards zero after a certain on- 
set temperature. The magnitude of the susceptibility 
in the constant region of the curve and the onset tem- 
perature of the decrease are both extremely sensitive to 
t he K/X ratio. Susceptibility curves reported in Refs. 
l7l l8l9lllCll for Eu'^+ compounds are very similar to those 
calculated within our model for small K. For larger val- 
ues of K, beginning with K/X ~ 0.1, the curve begins 
to develop a maximum at non-zero temperature which 
then evolves, with growing K, more toward a standard 
AFM susceptibility as shown in the inset of Fig. [Sjfor the 
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FIG. 9: The magnetic susceptibility of the two-site model 
with various levels of spin-spin interaction. The curve resem- 
bles a Van Vleck paramagnetic susceptibility for small K, and 
smoothly evolves to a shape of susceptibility typical of local 
AFM correlations when K dominates. 

limit K/X — !■ oo. Since the K/X ratio is manipulable by 
pressure, the effect of increasing K (signalling increasing 
spin-spin correlation) should be observable in the suscep- 
tibility curve, even for the conservative estimates of K/X 
mentioned earlier. 

The first principles calculation indicates that a non- 
zero spin-spin interaction of conventional strength does 
indeed exist. Though the energy difference between spin 
orientations for EuN is small at ambient volume, its mag- 
nitude can be increased with pressure. For EuP, the en- 
ergy difference is largest at equilibrium volume and is 
sensitive to manipulation by pressure. The model cal- 
culation shows that any non-zero interaction results in 
nonvanishing spin correlation {Si ■ Sj). 

VI. DISCUSSION 

We have initiated the study of the electronic structure 
of EuN and EuP specifically, and of the spin coupling and 
resulting behavior of the 4/^ ion more generally. The 
strong spin polarization of the 4/ shell, and the cou- 
pling of the moments L-|-S=J=0 to produce a nonmag- 
netic ion, makes a prediction of the physical electronic 
structure somewhat uncertain. We obtain, however, that 
both EuN and EuP are semimetallic due to overlap of 
the valence pnictide p bands with the conduction Eu 5c? 
bands. Recent measurement of the resistivity of pressed 
powders of EuN show a positive temperature coefficient, 
i.e. metallic or semimetallic conduction— 

We point out that even though the expectation value 
of any (vector) spin may be zero, exchange coupling 
will give rise to nonvanishing spin-spin correlations be- 
tween coupled ions. This coupling may lead to an un- 
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conventional phase transition in which the spin correla- 
tion length diverges but there is little or no signal in the 
magnetic properties. The determination of the order- 
ing temperature acquires new aspects: since entropy is 
(nearly) absent in both ordered and disordered states, a 
different mechanism from simple energy-entropy balance 
may control the ordering temperature. One possibility 
is that spin-phonon interactions promote disorder: rais- 
ing the temperature increases thermal vibrations, which 
modulates the interatomic separations and hence mod- 
ulates the exchange couplings dynamically. The phase 
transition would occur at a critical value of the exchange 
coupling "disorder." 

There are other Eu rocksalt compounds (divalent, 4/^) 
with FM ordering, as well as those with AFM ordering. 
In some cases, such as EuSe, a transition between the 
two magnetic alignments can be driven by manipulating 
the U parameter«ii or by applying pressure. Based on 
our LDA-t-U results, a state in which the spin moments 
of Eu ions are antiferromagnetically correlated is unlikely 
for EuN but probable for EuP. For EuN, a hidden ferro- 
magnetic correlation, or ordering, of spins (accompanied 
by an identical, but oppositely pointed arrangement of 
orbital moments) is energetically favored. It could be an 
experimental challenge to see the hidden "ferromagnetic 
order," which does not break any crystal symmetry. 
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1.55 
8.4 
1.11 
1.26 
0.41 
0.92 
0.37 



0.95 
1.11 
8.7 
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0.41 
0.48 



0.53 
1.26 
0.50 
9.4 
0.50 
1.26 
0.53 



0.48 
0.41 
1.90 
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1.11 
0.95 
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1.26 
1.11 
8.4 
1.55 



0.74 
0.37 
0.48 
0.53 
0.95 
1.56 
9.4 



IX. APPENDIX II 



VII. ACKNOWLEDGMENTS 

The authors are indebted to A. B. Shick for many dis- 
cussions on this work and for the code to produce the 
data presented in Appendix I. We acknowledge discus- 
sions with Z. Fisk, J. Lynn, S. M. Kauzlarich, J. Kunes, 
I. I. Mazin, D.A. Papaconstantopoulos, P. Novak, and D. 
J. Singh. W. E. P. was supported by DOE grant DE- 
FG03-01ER45876, DOE's Computational Materials Sci- 
ence Network, and the SSAAP program at LLNL under 
grant DE-FG03-03NA00071. 



VIII. APPENDIX I 

With U = 8 eV and = 1 eV as used in our LDA+U 
calculations, the direct and exchange matrix elements for 
diagonal occupation numbers are given by the following 
tables (the matrix elements are ordered from = -3 ^ iz 
— 3 from left to right and up to down): 



9.4 
7.8 
7.2 
7.1 
7.2 
7.8 
9.4 



7.8 
8.4 
7.9 
7.8 
7.9 
8.4 
7.8 



7.2 
7.9 
8.7 
8.4 
8.7 
7.9 
7.2 



7.1 
7.8 

8.4 
9.4 
8.4 
7.8 
7.1 



7.2 
7.9 
8.7 
8.4 
8.7 
7.9 
7.2 



7.8 
8.4 
7.9 
7.8 
7.9 
8.4 
7.8 



9.4 
7.8 
7.2 
7.1 
7.2 
7.8 
9.4 



We give here the parentage of the | J, Jz) states in terms 
of the the product states \S,Sz]L,Lz >= \Sz,Lz) {S = 
3,L — 3 are fixed. For our considerations only the J = 
and J — I states are needed. 



|0,0) 



|i,o) = 



1 

77 



(-3,3)-|-2,2) + |-l,l)-|0,0) 



|l,-l)-|2,-2) + |3,-3) 



a/3|2, -3) - V5|l, -2) + \/6|0, -1) 



2V7 



-V6\ - 1, 0) + \/5| - 2, 1) - %/3| - 3, 2) 



1 



^^^-3|-3,3)+2|-2,2)-|-l,l) 
+ |l,-l>-2|2,-2>+3|3,-3) 
^ ^\/3|3,-2) - \/5|2,-l) + V6|l,0) 



2V7 

-V6\0, 1) + \/5| - 1, 2) - V3| - 2, 3) 
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